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Abstract 

In this Letter we re-examine the idea that gravitational waves are required as a braking mechanism to 
explain the observed maximum spin-frequency of neutron stars. We show that for millisecond X-ray 
pulsars, the existence of spin equilibrium as set by the disk/magnetosphere interaction is sufficient to 
explain the observations. We show as well that no clear correlation exists between the neutron star 
magnetic field B and the X-ray outburst luminosity Lx when considering an enlarged sample size of 
millisecond X-ray pulsars. 

Subject headings: pulsars: individual(SAX J1808.4-3658, XTE J1814-338, ICR J00291+5934, XTE 
J1751-305) — gravitational waves 



1. INTRODUCTION 

Some accreting neutron stars in low mass X-ray bina- 
ries are thought to be spun up by transfer of angular 
momentum from an accretion disk. The spin-up pro- 
cess brings several neutron stars (NS) to spin in the mil- 
lisecond range. We refer to these systems as millisec- 
ond X-ray pulsars (MXPs), which comprise accretion- 
powered millisecond pulsars (with pulsations formed via 
magnetic channeled accretion), and nuclear-powered mil- 
lisecond pulsars (with burst-oscillations observed during 
thermonuclear bursts). These are very old systems with 
lifetimes of several billion years (van Paradijs & White 
1995). 

The spin-frequency {vs) of MXPs has been measured to 
date in 22 systems (Patruno 2010a, Papitto et al. 2011a) 
with a range between « 182 Hz and « 620 Hz. The 
narrow range of spin frequencies appears surprising con- 
sidering that the binary lifetime is billion years whereas 
the spin-up process operates on timescales of 10^ — 10* 
yr (White et al. 1988). 

Chakrabarty et al. (2003) and Chakrabarty (2008) ob- 
served that the spin-frequency distribution of these sys- 
tems (in 2003 comprising only 11 systems) was consistent 
with a uniform distribution with a cutoff frequency of 730 
Hz. Patruno (2010a) has shown how the same cutoff is 
found today with a sample size that is twice as large. The 
existence of a cutoff at such a relatively low frequency is 
somewhat unexpected because the Keplerian break-up 
frequency for a typical NS with any realistic equation of 
state is well above the 730 Hz cutoff. 

White & Zhang (1997) (henceforth WZ97) argued that 
the tight range of spin periods over two orders of magni- 
tude in observed luminosity, could be explained if these 
systems had reached spin equilibrium. In this picture 
the stellar magnetic field truncates the disk close to or 
at the corotation radius (where the stellar spin-frequency 
is equal to the Keplerian frequency of the disk), and no 
net angular momentum is transferred onto the star. 

This interpretation leads to an unexpected correlation 
between the NS magnetic field at equilibrium B and the 
outburst X-ray luminosity Lx- Indeed in transient sys- 
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tems Lx is given by the average mass accretion rate dur- 
ing an outburst, which depends on the outburst length 
and on the outburst recurrence time. These quantities 
are set solely by the accretion disk physics (Lasota 2001), 
so it seems difficult to justify a correlation with the mag- 
netic field. 

Following an earlier suggestion of Papaloizou & Pringle 
(1978) and Wagoner (1984), Bildstcn (1998) proposed 
a different mechanism to limit the spin-up of MXPs, 
introducing a loss of angular momentum via gravita- 
tional wave (GW) emission. In this scenario a corre- 
lation between B and Lx is much weaker and a sharp 
spin-frequency cutoff is expected due to the strong spin- 
frequency dependence (oc v^) of the spin-down torque. 

However, recent observational results on some MXPs 
have shown that the efficiency of GW induced spin-down 
for MXPs might be very low. Patruno (2010a), Hart- 
man et al. (2011) and Papitto et al. (2011b) measured 
the long-term evolution of the spin-frequency of IGR 
J00291-I-5934, revealing both accretion torques during 
one outburst and a long term spin-down during quies- 
cence, which they interpreted as due to magneto-dipolc 
torques. The combined effect of spin- up in outburst 
and spin-down in quiescence leads to a timescale of sev- 
eral billion years for the evolution of the 599 Hz spin- 
frequency of this MXP (Patruno 2010a). Hartman et al. 
(2011) and Papitto et al. (2011b) also showed that GW 
torques are unimportant for MXPs with frequencies up 
to 599 Hz. 

Another recent result (Hartman et al. 2008, 2009, 
HaskeU & Patruno 2011), showed that in MXP SAX 
J1808.4-3658 and XTE J1814-338 stringent upper hmits 
can be set on the accretion- induced spin- up of the NS. 
Haskell & Patruno (2011) in particular, rejected the sce- 
nario in which a continuous train of GWs is responsible 
for balancing the spin-up process. Those authors pro- 
posed spin equilibrium as set by the disk/magnetosphere 
interaction as the most likely mechanism to explain the 
observed behavior of the two MXPs. 

In this Letter we revisit the analysis of WZ97 and show 
that: 

• the observations suggest that a large majority of 
MXPs can be in spin equilibrium without implying 
a correlation between the B field and the X-ray 
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outburst luminosity 

• invoking GWs to explain the cutoff in the spin dis- 
tribution is not necessary given the uncertainty in 
models of more realistic disk/magnetosphere inter- 
actions 

2. MILLISECOND X-RAY PULSAR SAMPLE 

We use 18 MXPs (sec Tabic 1) which have both con- 
firmed spin frequencies (see Patruno 2010b and refer- 
ences therein) and a realistic determination of the out- 
burst flux (see Watts et al. 2008, Riggio et al. 2011). 

We do not use any accreting NS whose Vg is inferred 
only from observations of twin-kHz QPOs. The idea 
that the separation A v between twin-kHz QPOs is di- 
rectly related with Vs was supported by the observed 
ratio /S.v/vs ~ 0.5 — 1.0 in a few MXPs (van der Klis 
et al. 1996, Wijnands et al. 2003). It has recently been 
shown, however, that twin-kHz QPOs might not be a 
good proxy for NS spin frequencies and could even be 
completely unrelated to Vg (Mendez & Belloni 2007, Yin 
et al. 2007). 

In comparison, WZ97 used a sample of 10 accreting 
NS, five of which had spin frequencies inferred only from 
the observation of twin-kHz QPOs and five from burst- 
oscillations. 

3. DISK/MAGNETOSPHERE INTERACTION 

The interaction between a star's magnetic field and 
the surrounding accretion disk provides a natural mech- 
anism to limit the spin-frequency of the star. By assum- 
ing that most systems have reached their limiting spin- 
frequency, WZ97 estimated the magnetic field strength 
for the sources in their sample and found an unexpected 
correlation between B and Lx ■ 

For strong magnetic fields, the disk will be truncated 
at some distance from the star (at the magneto spheric 
radius, rm), which is generally estimated to be (Pringle 
& Rees 1972): 

r™ - i^iMZlC M^l'' RfJ' bI'' km (1) 

where M_io is the accretion rate in units of lO~^''A/0/yr, 
Ml. 4 is the NS mass in units of 1.4 Mq, i?io its radius 
in units of 10 km and ^ parametrizes the uncertainties 
in evaluating the torque at the edge of the accretion 
disk (thought to be in the range ^ w 0.3 — 1, Psaltis 
& Chakrabarty 1999). 

Inside rm gas will be channeled onto the star, spinning 
the star up at a rate: 

j = M^GMrra. (2) 

However, if the star is spinning very fast, the disk may 
be truncated outside the corotation radius r^, where the 
Keplerian frequency is equal to the spin-rate of the star: 

rco = lQ^'iMl'2vs'^'^kuY, (3) 

and the spinning magnctosphere will present a centrifu- 
gal barrier that can inhibit accretion (the so-called "pro- 
peller regime", Illarionov & Sunyaev 1975). This sug- 
gests that the star will eventually spin at a rate such 
that rm — Tco ( "spin equilibrium" ) . 



Under this assumption, the magnetic field at equilib- 
rium can be estimated as: 

i? = 8.8 X IQi" r^'^M^Jw^ifiR^^y:'"' G. (4) 

The ten accreting NS used by WZ97 showed a spread 
in luminosities over two orders of magnitude whereas the 
spin periods clustered between 2.8 and 3.8 s (263—362 
Hz). Equation 4 shows that the small range of spin pe- 
riods over a large span in luminosity can be explained 

1/2 

within the spin equilibrium scenario only if i? oc L-^ 

(under the reasonable assumption that Lx cx M) . How- 
ever, the sample used by WZ97 suffered of two serious 
biases: 

• all bright systems have NS spin frequencies deter- 
mined only from twin-kHz QPOs, which today is 
not considered a robust method 

• in two cases the spin period was taken at twice the 
observed value to match the /S.v of twin-kHz QPOs. 

If we remove the first bias by ignoring twin-kHz QPO 
sources, then the luminosities span less than one order of 
magnitude instead of two. By removing also the second 
bias by using the observed spin period (instead of twice 
its value), the spin periods span a range between 1.7 and 
2.8 ms (362—589 Hz). The difference between the min- 
imum and maximum spin frequency therefore becomes 
more than twice the value used by WZ97. Therefore 
there is no clustering of MXP spin frequencies over two 
orders of magnitude in luminosity as reported by WZ97. 

We can go further and repeat the analysis of WZ97 
on the enlarged sample of MXPs given in Table 1, with 
the spin periods spanning now an even broader range 
between 1.6 and 5.5 ms (182-620 Hz). Assuming each 
source is in equilibrium and the magnetic field is given 

1/2 

by Eq. 4, we see no convincing correlation B oc L-^ 
(see Figure 1). The field at equilibrium has now a large 
spread which depends on the broad scatter in observed 
Vs- Fig. 1 shows also that the field is not stronger for 
bright sources and weaker for faint sources as was inferred 
by WZ97 (compare with Figure 1 in WZ97). 

4. REALISTIC MODELS OF DISK/MAGNETOSPHERE 
INTERACTIONS 

The calculation in § 3 does not incorporate the consid- 
erable uncertainties in calculating how angular momen- 
tum can be extracted from the star, which will strongly 
affect the spin equilibrium frequency. Here we revisit a 
more realistic model for the interaction between an ac- 
cretion disk and the magnetic field, and demonstrate how 
these uncertainties can lead to a large variation in equi- 
librium spin periods for a given B and M . 

The estimate for the magnetosphcric radius given by 
Eq. 1 does not consider the relative rotation between the 
star and the disk, and as such is only applicable when 
the disk is truncated well inside the corotation radius 
(that is, the magnetic field is rotating relatively slowly 
compared with the inner parts of the disk). As the disk 
moves outward, the relative rotation between the disk 
and star becomes critical in determining the location of 
r„i and the disk structure. Once r,„ is very close to 
Tco, Eq. 1 fails completely, since now the magnetic field 
presents a centrifugal barrier to accretion. D'Angelo & 
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TABLE 1 

Millisecond X-ray Pulsar Sample With Spin-Frequency and Luminosity 



Source Name 


Lx 


Spin Frcq. 


Outb. 


Flux 


Distance 




(I0^%rg/s) 


(Hz) 


(10-" 


erg/s/cm"^) 


(kpc) 


4U 1608-522 


40 


620 


2.0 




4.1(0.4) 


Aql X-1 


32 


550 


1.3 




4.55(1.35) 


SAX J1748. 9-2021 


31 


442 


0.4 




8.1(1.3) 


KS 1731-260 


30 


524 


0.49 




7.2(1) 


XTE J1751-305 


28 


435 


0.29 




9(3) 


IGR J17191-2821 


18 


294 


0.26 




7.5(3.5) 


IGR J17511-3057 


15 


244 


0.2 




8(2) 


MXB 1659-298 


13 


567 


0.08 




12(3) 


XTE J0929-314 


7 


185 


0.1 




7.8(4.2) 


SAX J1750. 8-2900 


7 


601 


0.12 




6.79(0.14) 


GRS 1741.9-2853 


6 


589 


0.1 




7.2(2.8) 


XTE J1807-294 


6 


191 


0.072 




8.35(3.65) 


SAX J1808. 4-3658 


5 


401 


0.35 




3.5(0.1) 


IGR J00291+5934 


5 


599 


0.16 




5(1) 


XTE J1814-338 


4 


314 


0.069 




6.7(2.9) 


SWIFT J1756. 9-2508 


3 


182 


0.04 




8(4) 


HETE J1900. 1-2455 


2 


377 


0.09 




4.7(0.6) 


NGC6440 X-2 


2 


206 


0.019 




8.1(1.3) 



WZ97 model 
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Fig. 1. — Left Panel: Inferred magnetic field at equilibrium vs. luminosity relation for MXPs as reported in Table 1. The black 
dots refer to magnetic fields inferred under the hypothesis of spin equilibrium Vm = Tco- The Eddington luminosity for a 1.4 Mq is 
LecLcL = 1-8 X 10^*^ergs~^. Right Panel: Inferred magnetic field at equilibrium vs. spin period. The plot uses the same sources and 
symbols as in the left panel. The three lines refer to luminosities of 0.1%, 1% and 10% Eddington. 



Spruit (2010) demonstrated that for ~ Tco, the an- 
gular momentum added by the disk-field interaction can 
change the disk structure and lead to a disk truncated 
just outside r^o without any accretion onto the star or 
outflow. 

Calculating the spin rate in equilibrium for a given ac- 
cretion rate is thus not as simple as equating Eq. 1 with 
Eq. 3, and requires a good understanding of the spin- 
down torques, which are subject to large uncertainties. 
For the picture presented in D'Angelo & Spruit (2010) 
the spin-down torque comes from the interaction between 
the disk and field outside co-rotation. This leads to a 
spin-down torque of magnitude: 



where ^ = B^R^ is the magnetic dipole moment, rj ^ 1 
is a dimensionless coefficient that quantifies the strength 
of the disk-field coupling, and ^r/rm < 0.3 indicates the 
radial extent of the coupling between the star and the 



disk. Both these parameters will be dependent on the 
detailed disk-field coupling, and may vary between sys- 
tems as a result of, e.g., the magnetic field inclination 
or the large-scale magnetic field configuration (D'Angelo 
& Spruit 2011). An additional complication in predict- 
ing the amount of spin-down torque is the appearance 
of episodic cycles of accretion (Spruit & Taam 1993; 
D'Angelo & Spruit 2010, 2011), in which mass accumu- 
lates just outside rc and is periodically dumped onto the 
star. D'Angelo & Spruit (2011) have found that the net 
amount of spin-down in cyclic accretion can be up to 
100% larger than steady accretion at the same accretion 
rate. 

Other sources of angular momentum loss from the star 
can come from magneto-dipolc radiation or stellar out- 
flows. Magneto-dipole radiation may dominate the spin- 
down for sources in quiescence (Hartman et al. 2008, 
2009), but is unlikely to dominate over the disk-field 
interaction while the star is in outburst. In contrast, 
a considerable amount of angular momentum could be 
expelled via a stellar wind, depending on how much 



4 




Fig. 2. — Allowed parameter space (shaded gray area) between 
the inferred magnetic field at equilibrium and the source luminos- 
ity for a hypothetical 730 Hz accreting NS. The top and bottom 
dashed lines are obtained by equating the spin-down in Eq. 5 and 
the spin-up in Eq. 2 and are calculated for r?^r- = 10"'^ and 0.3, 
respectively. The solid lines refer to the radiation pressure domi- 
nated disk solution of Andersson et al. (2005), with ^ = 0.3 and 
1. For bright systems the inferred B field at equilibrium has an 
uncertainty that spans almost three orders of magnitude. 

mass is ejected in an outflow (Matt & Pudritz 2005). 
Stellar winds have been observed in numerical simula- 
tions to carry away a large amount of angular momen- 
tum, and could compete or even dominate over field-disk 
coupling as a source of angular momentum loss in the 
star(Romanova et al. 2009). 

Finally, at high accretion rates, the coupling between 
the disk and field can also change the amount of angular 
momentum deposited onto the star from accretion (An- 
dersson et al. 2005). In this case the inner parts of the 
disk become radiation-pressure dominated, and the spin 
equilibrium condition translates then into the relation: 

Bo^fl-^) lTi^:^^ (6) 

V i^EddJ 

(see Eq. 38 and 45 in Andersson et al. 2005). 

Summarizing, the value of jy-^ can be different in dif- 
ferent systems and this might introduce a spread of B 
fields at equilibrium up to an order of magnitude. For 
bright sources this spread can increase up to three orders 
of magnitude because of the effect of radiation-pressure 
dominated disks. Therefore we would not expect to ob- 

1/2 

serve B oc L-^ even for sources that were all spinning 
at the same frequency. We show this in Fig. 2 where 
we have chosen an hypothetical accreting NS spinning 
at 730 Hz and a range for the parameter 77— in Eq. 5 
between 0.001 and 0.3 (D'Angelo & Spruit 2011). We 
have also included the possibility that the disk becomes 
radiation pressure dominated as described in Andersson 
et al. (2005), with £, = 0.3 and 1. 

The allowed region in Fig. 2 shows that a B field of 
~ 10^ G is already sufficient for an MXP spinning at 
730 Hz at equilibrium regardless of the source luminos- 
ity. Furthermore, since B (x i^s any accreting NS 
spinning above 730 Hz also requires fields of the order of 
10* G or lower at equilibrium. Therefore there is no need 



for bright systems like Sco X-1 to have larger B fields at 
equilibrium than faint sources. 

5. DISCUSSION 

In this Letter we have re-assessed the possibility that 
the apparent cutoff at 730 Hz in the spin distribu- 
tion of accreting NS may be due to these systems 
being at or close to spin equilibrium, as set by the 
disk/magnetosphere interaction. This possibility was ex- 
amined in detail by WZ97, who concluded that such 
a scenario would require an unexpected correlation be- 
tween the magnetic field of the star and the outburst lu- 
minosity. Their analysis was, however, biased by the rela- 
tive scarcity of observations and by the inclusion in their 
sample of bright systems for which the spin-frequency 
was inferred from the separation of the twin-kHz QPOs, 
a method which is now known to be unreliable (Mendez 
& Belloni 2007). Nevertheless, the work of WZ97 did lead 
to the suggestion that the disk/magnetosphere interac- 
tion alone was insufficient to explain the spin-frequency 
distribution of accreting NS, unless the B field at equi- 
librium was correlated with the X-ray luminosity of the 
binaries. To avoid this correlation, Bildsten (1998) pro- 
posed that GWs may remove angular momentum at a 
rate sufficient to explain the observed spin distribution. 

Repeating the analysis of WZ97 with the extended 
sample of systems available today and removing the twin- 
kHz QPO systems, we have shown that no strong corre- 
lation between the magnetic field strength and luminos- 
ity is required to explain the observed spin distribution. 
This is due to two main effects. First the spread in ob- 
served spin frequencies has considerably increased since 
the work of WZ97, leading to a larger scatter in the val- 
ues for the magnetic field strength required for equilib- 
rium, even in the simplest disk/magnetosphere picture. 
On the other hand the disk structure itself and the source 
of spin-down torque vary from system to system, leading 
to a considerable range of equilibrium B for a given 1/5 
and Lx- 

From an observational point of view, there is evidence 
that two of the MXPs may be close to spin equilibrium 
during the outburst (Haskell & Patruno 2011), while 
the MXP IGR J0029H-5934 (and to some extent XTE 
J1751-395, Riggio et al. 2011) exhibits a saw-tooth like 
behaviour, with the spin-up during the outburst nearly 
balanced by the spin-down during quiescence (Patruno 
2010a Papitto et al. 2011b, Hartman et al. 2011). This 
behaviour suggests these systems arc in long-term spin 
equilibrium, (Eisner et al. 1980, Lamb & Yu 2005), a sce- 
nario which would be consistent with the inferred mag- 
netic field of IGR J00291-h5934 (B w 2 x 10* G) and 
XTE J1751-305 (B « 4 x 10* G). 

Furthermore, Haskell & Patruno (2011) have shown 
that GW emission is not consistent with the timing prop- 
erties and quiescent luminosity of SAX J1808. 4-3658 and 
XTE J1814-338. Note that Ho et al. (2011) have sug- 
gested that GW emission may be required to explain the 
lack of observed systems with low spin rates and long 
orbital periods. This argument is, however, based on the 
assumption that the spin-up torque is of the simple form 
in § 3 and might be necessary to reassess it with more 
realistic models. 

The situation may be quite different for bright persis- 
tent sources like Sco X-1 and other persistent GX and 
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Z sources (see van der Klis 2004 for a review on the 
source classification). For most of these sources the spin- 
frequency of the NS is not known, and the high accretion 
rate would imply spin equilibrium, for the simple accre- 
tion torque of § 3, well above the Kcplcrian breakup 
frequency for magnetic field strengths of the order of 
B « 10* G. Although it is possible that these systems 
do contain a sub-millisecond pulsar, this might be dif- 
ficult to reconcile with the fact that the distribution of 
millisecond radio pulsars also appears to have a cutoff 
at around 700 Hz (Hessels et al. 2006). However, unlike 
X-ray pulsars, the sensitivity of radio surveys to sub-ms 
pulsars degrades with the unknown dispersion measure, 
which might partially explain why radio sub-ms pulsars 
have never been found so far (Hessels et al. 2007). 

A second possibility is that these systems have reached 
spin equilibrium at ^ 700 Hz with a stronger magnetic 
field in bright systems, as predicted by Konar & Bhat- 
tacharya (1999). If the B field is confined to the NS 
crust a positive correlation exists between the average 
mass accretion rate (which, unlike the outburst accre- 
tion rate, does not depend on the details of the disk 
instability model) and the final field strength. In this 
scenario, the mechanism responsible for the B field de- 
cay is Ohmic dissipation in the NS crust, which is sub- 
stantially accelerated with accretion. The accreted mate- 
rial also pushes the current producing the B field deeper 
down to the crust-corc interface, where resistivity is low 
and Ohmic dissipation becomes unimportant (see e.g.. 



Romani 1990). 

An alternative scenario is that these bright systems 
may, in fact, be emitting GWs. The higher tempera- 
tures of these systems and the constant stream of ac- 
creted material make it likely that the star could sustain 
a mountain in the crust, due to compositional and heat- 
ing asymmetries (Ushomirsky et al. 2000) or that the 
magnetic field may be distorted by the accretion flow 
and support a substantial quadrupole (Melatos & Payne 
2005). 

However, a plausible minimal scenario is that, as dis- 
cussed in § 4, the spin-up torque is much weaker than the 
estimate in § 3 at higher luminosities, while the magnetic 
field strength is still in the region of i? « 10* G. This 
could either be due to the inner parts of the disk being ra- 
diation pressure dominated (Andcrsson et al. 2005), the 
disk/magnetosphere coupling being weaker (D'Angelo & 
Spruit 2010) or to additional sources of spin-down, such 
as mass outflows (Romanova et al. 2009). 

It would thus appear that the observed spin distribu- 
tion is consistent with the notion that all systems are 
either at, or close to, spin equilibrium. This picture does 
not require a strong correlation between the magnetic 
field strength and the luminosity and is compatible with 
modest field strengths of the order 10* G, with no need 
to assume much higher field strengths for any systems. 
This field strength is compatible with the B field that is 
inferred from the quiescent spin-down of three MXPs. 
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